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Abstract 

We reconsider the evolution of strongly degenerate neutrinos in the early universe. Our chief 
concern is the validity of the entropy conservation after the neutrino annihilation process has 
frozen out (so that the establishment of chemical equilibrium is not trivial). We argue that the 
entropy indeed conserves because elastic scattering keeps the neutrino and antineutrino distribution 
functions in the equilibrium form and the sum of their chemical potential keeps zero even after 
the neutrino annihilation freeze-out. We also simulate the evolution of the degenerate neutrino 
spectrum to support the argument. We conclude that the change in the neutrino degeneracy 
parameter when the relativistic degrees of freedom in the universe decreases is calculated using 
the entropy conservation and the lepton number conservation without worrying about at what 
temperature the neutrino annihilation process freezes out. 
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I. INTRODUCTION 



The cosmology with the strongly degenerate neutrinos or, in other word, large lepton 
asymmetry is one of well-investigated themes. There are researches concerning how to 
generate it and what kind of cosmological consequences it produces at various cosmological 
epochs (see Ref. [lj). In this paper, we make comments on thermodynamic properties of 
such degenerate neutrinos in the early universe. Some of them have been already pointed 
out in Refs. and we make complementary arguments. Others provide corrections to 

probable misunderstandings found in Ref. Q]. 

First of all, we would like to compare Refs. and this paper. In Ref. it is argued 
that the evolution of the neutrino degeneracy parameter is given by imposing the neutrino 
number conservation after the neutrino annihilation freeze-out. Therefore the freeze-out 
temperature and the variation of the temperature are regarded to be necessary, where the 
latter is calculated by integrating the covariant energy conservation law because the entropy 
conservation is not considered to hold. As for the antineutrino degeneracy parameter £ p , 
the relation + £ p = is considered to break down. In Ref. Q], contrary to Ref. it 
is argued that the temperature variation can be calculated using the entropy conservation 
and that the relation + £p = holds after the freeze-out. There, the entropy conservation 
is regarded to hold approximately and the justification for its application is attributed to 
the exponentially suppressed antineutrino number density under the existence of the strong 
degeneracy. In this paper, we note that the entropy conservation is derived using the relation 
£v + £p = 0, which holds also after the freeze-out as described in Ref. Q]. Then the entropy 
conservation is valid at any time. Together with the lepton number conservation, we can 
calculate how the degeneracy evolves without knowing at what temperature the freeze-out 
takes place. This is the point we have missed in the previous paper and want to stress in this 
paper. However, we note that the calculations in Ref. |2| are not wrong because the neutrino 
number conservation used there is virtually same as the lepton number conservation when 
the degeneracy is large. Just we have done redundant calculations on freeze-out temperature. 

In Sec. m we first describe how the variation in the neutrino degeneracy parameter is 
calculated with the lepton number conservation assuming the entropy conservation. Then 
we show the entropy to conserve as long as neutrinos are kinematically coupled to the rest 
of the cosmic plasma. In Sec. IIII[ we show the results obtained by numerically simulating 
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the degenerate neutrino spectrum, 
equilibrium and £ v + — holds, 
conclusion. 



In this way, we directly confirm they are in thermal 
In Sec. IIVI we summarize the discussion and give the 



II. THE NEUTRINO DEGENERACY VARIATION CALCULATED BY THE EN- 
TROPY CONSERVATION 

The lepton asymmetry in the universe is of course measured by the lepton number density: 
ni = n u — rip. When neutrinos and antineutrinos are in thermal equilibrium, their number 
densities are specified by temperature T together with chemical potential p (or as is fre- 
quently used, with degeneracy parameter £ = fi/T) and annihilation process u+v <-> e _ + e + 
ensures (, u = — £ p so that ni = £{(£/7r) 2 + l}T 3 /6 (we assume for simplicity £ u > and denote 
it as £). 

Since the lepton number jVj = a 3 ni and the total entropy S = a 3 s in the universe conserve 
(where a is the scale factor and s is the entropy density), it is useful for quantifying the 
lepton asymmetry to consider their ratio rji = ni/s, which also conserves. Using energy 
density p, pressure P and number density n, s is calculated as s = (p + P — fin)/T, and is 
often written as s = (2tt 2 / '45) gs(£, T)T 3 where g s denotes the relativistic degrees of freedom. 
Then rji is related to the degeneracy parameter as, rji oc (£ 3 + 7r 2 £)/ 9s(£,,T). 

It is obvious from the last expression that £ takes different values as the total relativistic 
degree of freedom of the universe changes and how much it does can be calculated with the 
r]i conservation (especially, £ stays constant while total gs does not change and this is the 
reason that £ is often used to quantify the degeneracy). 

However, when the degeneracy is very large, there are occasions that the entropy con- 
servation is not trivial as explained below (the lepton number conservation evidently holds 
because it is respected in every relevant elementary process). With neutrino degeneracy, 
there exists more neutrinos and less antineutrinos than without it. Since this makes harder 
for neutrinos to find partners of annihilation process, it freezes out (the process rate becomes 
less than the cosmic expansion rate) at higher temperature. When the degeneracy is so large 
that neutrino annihilation freezes out before the muon and antimuon annihilate, we can not 
in general expect the entropy to conserve during the muon annihilation because the process 
which ensures the chemical equilibrium to hold has frozen out. 
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The higher neutrino annihilation freeze-out temperature under the existence of degener- 



acy has 



xhmi noticed in Rot". Q| and its effect on cosmology has lxvn discussed there and also 



in Ref. jj. However, as pointed out by Ref. they have regarded the neutrino annihila- 
tion freeze-out as the neutrino decoupling so they have concluded that neutrinos would not 
be heated while the muon annihilation. This is not true because although neutrinos come 
short of annihilation partners, antineutrinos, they have enough elastic scattering partners, 
for example electrons. The correct picture is that even the neutrino annihilation freezes out, 
they keep contact with the rest of the cosmic plasma and decoupling does not take place. 
In other words, after the neutrino annihilation freeze-out, the chemical equilibrium in not 
ensured to hold but the kinetic equilibrium is. 

In this picture, what occurs during muon annihilation is that neutrinos preserve number 
and keep equilibrium distribution with the temperature following photon's which falls more 
slowly than a -1 . That is, 

N v = a 3 n„ oc a 3 T 3 f - x = const. (1) 

J exp(y - + 1 

while aT deviates from unity and increases. Therefore, as is easily seen, £ has to decrease. 
On calculating how much aT increases and £ decreases, Ref. jl| has suggested the use of 
covariant energy conservation law because they consider the entropy conservation does not 
hold when there is degeneracy. Also, they have stated that the relation + £ p = 0, which is 
true during chemical equilibrium, breaks down. In summary, from Ref. [l|'s point of view, 
the evolution of £ would be calculated as: 1) calculate the neutrino annihilation freeze- 
out temperature, 2)use the total entropy conservation before the freeze-out, and 3)use the 
neutrino number conservation and the covariant energy conservation law after the freeze-out. 

In our recent paper, we expressed different opinions from theirs. In Ref. Q], we argued 
that the total entropy conserves and £„ + = holds even after the freeze-out. Our 
argument is as follows. 

The variation of the entropy is determined by the second law of thermodynamics j^j]: 

TdS = d(pV) + PdV - fidN. (2) 

The first two terms on the right-hand side vanishes according to the covariant energy con- 
servation. For the last term, there are contributions from neutrinos and antineutrinos so 
that 

TdS = -fi u dN u - fipdNp = -(fi u + fi p )dN u , (3) 



where for the second equality, we use the lepton number conservation dNi = d(N u — N p ) = 0. 
Before the neutrino annihilation freezes out, the particles are in chemical equilibrium so 
f^u + \i>v = holds and the entropy conserves even if dN v ^ (during for example muon 
annihilation). 

After the neutrino annihilation process has frozen out, since the chemical equilibrium 
breaks in general, we can not apply \x v + fip = immediately so it is not trivial that the 
entropy conservation holds when dN v ^ 0. However, what we call here "freeze-out" is for 
neutrinos and not for antineutrinos who have a lot of annihilation partners. As a result, the 
chemical potential of antineutrinos is expected to keep the value \i v = —fi v and the entropy 
conserves. 

To make this argument clearer, we consider the Boltzmann equation for the annihilation of 
antineutrinos during the muon annihilation. For that purpose, it is important to notice that 
even after the neutrino annihilation freeze-out, the elastic scattering is sufficiently frequent 
so every particle species are in kinetic equilibrium with certain well-defined temperature and 
as the universe expands, the kinetic equilibrium is maintained as the temperature decreases. 
The relevant part of the Boltzmann equation is 

dn v f d 3 p 1 d 3 p 2 d 3 p 3 d 3 p 4 4 (4) 

-dT + 3Hn > = J 2^(2^)3 2^(2^)3 2^(2^)3 2^(2^)3 (2?r) * ^ + P * ~ P * 

x|A^| 2 {[l - /p(£i)][l - f„(E 2 )]f e -(E 3 )f e+ (E 4 ) 

-HE^UE^l - f e -(E 3 )][l - / e +(£ 4 )]}, (4) 

where \A4\ 2 is the (angular integrated) invariant amplitude squared. Suppose that at first 
the universe has temperature T\ with \i v + \i v = 0. At that moment, the number density 
obeys the Boltzmann equation without the expansion term dn D /dt = 0. This is consistent 
with the collision term, the right hand side of Eq. (@J, which vanishes when the distribution 
functions take equilibrium form f(E) = l/[exp{(E — /i)/Ti} + 1] with /i e ± = 0, \i v + fip = 
and energy conservation E\ + E 2 = E3 + E±. As the muons annihilate and the universe 
expands, the equilibrium with temperature Ti breaks but by the frequent elastic scattering, 
it quickly settles down to next equilibrium with T 2 (< Ti). This transition is driven by the 
Boltzmann equation like Eq. (j3J) with electron and positron temperature T 2 and neutrino 
and antineutrino temperature T\. The difference takes place because the electromagnetic 
interaction is much stronger than the weak interaction and electrons are immediately heated 
by annihilating muons but neutrinos are not. This gives non-zero collision term to evolve the 
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distribution functions or, knowing they take equilibrium form, to evolve chemical potentials. 
The relation between the chemical potentials is determined by demanding the Boltzmann 
equation be drip/dt = when H — > and T% — > T 2 . This is only achieved when the phase 
space factors on the right hand side of Eq. (0J cancel out i.e. /!„ + ii v — 0. 

We would like to note three points. First, similar argument does not work with neutrino 
Boltzmann equation because after the neutrino annihilation freeze-out, when collision term is 
negligible to expansion term, it is dn u /dt + 3Hn u w and does not give us information about 
distribution functions. Second, the argument above assumes the period in which temperature 
decrease is dictated by the muon annihilation in addition to the cosmic expansion. On the 
other hand, when particle degree of freedom is constant, the universe becomes cooler as 
a whole only by the cosmic expansion so there occurs no temperature difference between 
electrons and neutrinos. Then the terms in the right hand side of Eq. (j3J) cancel out to 
give the evolution equation dn D /dt + 3Hn p properly expressing the number conservation. 
Third, it is crucial that there exists temperature difference between electrons and degenerate 
neutrinos to show \i v + fip = but the difference is quickly erased due to the frequent elastic 
scattering and we can regard T v = T 7 for the cosmological time scale. More concretely, 
the argument in the previous paragraph is valid when T\ — T 2 is much smaller than the 
temperature difference between before and after the muon annihilation. 

Now that we show \i v + /i p = 0, it is readily seen from Eq. Q that the total entropy 
conserves even after the neutrino annihilation process freezes out. Then, the evolution of £ u 
is calculated with the total entropy conservation and the lepton number conservation with 
£p = — £„. Note that the neutrino annihilation freeze-out temperature is not necessary for 
the calculation, contrary to what is discussed in the literatures such as Refs. jlj] 0] . The 
naive treatment we have introduced at the beginning of this section turns out to be correct. 

III. NUMERICAL SIMULATION OF THE NEUTRINO SPECTRUM EVOLTU- 
ION 

In this section, we simulate the evolution of degenerate neutrino spectrum to confirm 
the argument given in the previous section. We find manifestly the thermal equilibrium 
distribution is preserved with the same temperature as the photons and the relation £ v +£& = 
holds . 
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Before showing the results, we describe our simulation method. Similar simulations are 
performed in Refs. 

BB 

and more details are found. We assume for simplicity that only 
electron-type neutrinos are degenerate and other types have no degeneracy. We evolve neu- 
trino and antineutrino distribution functions, f u (y) and fp(y), and the photon temperauture 
T as functions of x = m Q a where m is an arbitrary energy scale (we use m = 1 GeV and 
set it unity hereafter), y is defined by y = ap where p is the particle momentum. 
With these variables, derivatives of the distribution functions are calculated by 

dfv{v) _ Cv{v) /_\ 

dx ~ ~Hx~ : 1 ' 

where C is the collision term and H is the cosmic expansion rate. H is calculated from the 
total energy density ptot as H = y/ptot/3/Mpt where M pl = 2.436 x 10 18 GeV is the Planck 
energy. 

For the collision term, we include the elastic scattering v{y) + <-> v{y) + e ± and the 
annihilation v + v <-> e" + e + . We denote the former contribution as C e and the latter C a 
so that C = C e + C a . We calculate them with approximation for electrons to be massless 
and to obey Boltzmann statistics as in Ref. Then 
2G 2 F [(C V + l) 2 + (C A + l) 2 ]x 



n 3y2 



dy'[l-f v (y')}F l (y,y') + J dy'[l - f u (y')}F 2 (y,y' 

( ry r°° ^ 

+[l-f v (y)][j o dy'W)B 1 (y,y') + J dy'f v (y')B 2 ( y , y ') j j , (6) 

where Gp is the Fermi coupling constant, Cy = —1/2 and Ca = — 1/2 + 2sin 2 # (6: weak 
mixing angle). Functions F and B are defined by 



Fi{y,y') 


= D(y,y') + E(y,y')exp 




= D(y,y)e W ( ^ ) 


Bi(y,y') 


= F 2 (y\y), 


B 2 {y,y') 


= Fi{y',y), 



y_ 

xT 



(7) 



E(y,y')[-i), (S) 



xT 



(9) 
(10) 



where 



D(y,y') = ±- r {y* + y >* + 2(y- y ')xT + 4x 2 T 2 }, (11) 
E(y, y') = -—^ {y 2 y' 2 + 2yy'(y + y')xT + 2(y + y'fx 2 T 2 + A(y + y')x 3 T 3 + 8x 4 T 4 }(12) 



As for C a , the expressions in Ref |8( have not included neutrino degeneracy so we need 
to modify the annihilation term to be 

4G|[(CV + l) 2 + (C A + l) 2 ] 



y + y' 



97T 3 X E 

x I dy'yy' 3 [U{v)Ml/) ~ [1 - UvW ~ Ml/)] e*P ) "( 13 ) 

For antineutrinos collision terms, we have to just exchange f v and f p in Eqs © and (|13j). 

The derivative of the temperature is obtained from the covariant energy conservation 
dp/dx = —3(p + P)/x. Since our simulation includes photons, electrons (approximated to be 
massless), muons, two types of neutrinos with no degeneracy and one type with degeneracy, 

d T _ H 9x 4 P7 + sfo, + p,)} + I d yy"t + 2^1 



dx 



dx ~ 9 x 3et + 3fi£i 

dT ^ dT 



(14) 



As the initial condition, we take a; = 10 which corresponds to T = 100 MeV and assume 
degenerate neutrinos have equilibrium distribution with £„ = — £ p = 10. Actually, for this 
condition, since the annihilation rate exceeds cosmic expansion rate (see Ref. Q), nentrinos 
are in chemical equilibrium. The muon mass is about 106 MeV so they are almost fully 
relativistic (relativistic degree of freedom is 3.11, while fully relativistic particle would have 
3.5) at 100 MeV. We follow the evolution down to until about 10 MeV at which temperature 
the muons almost annihilate away. 

We note here some technical detail concerning numerical calculation. To discretize mo- 
mentum, we take equally spaced 100 points in < y < 20. Time step is fixed to Ax = 10~ 4 . 
Since the differential equation for f's are stiff but not for T, /'s are first evolved with 
2nd order semi-implicit method and then, using that results, T is evolved with 2nd order 
Runge-Kutta method. 

The simulation results are summarized in Figs. ^ to H3 Figs. ^ and ^indicate the spectra 
keep equilibrium form with well-defined temperature and chemical potential. Fig El shows 
how the degeneracy parameter evolves which can be reproduced using the entropy conser- 
vation and the lepton number conservation. 

For Figs. H] and what we actually evolve is the distribution functions f v (y), fp(y) 
and photon temperature T but we express the results in terms of momentum dependent 
degeneracy parameters £„(y) and £p(y) which are calculated using the relations 

1 1 
exp( — + 1 ex Pl^r t,v{v)) + 1 
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FIG. 1: This figure shows the evolution of £ v (y). £ is computed as in the text. The appearance 
of the horizontal lines indicates that the neutrino distribution takes equilibrium form at each 
temperature. We plot £ v at several temperatures. The greater line intervals correspond to the 
faster decrease in the muon relativistic degree of freedom. 

or €u(p) — y/( x T) — ln(/~^ — 1). Note that we do not assume neutrinos and antineutrinos 
have equilibrium distribution at this stage. The result that £'s computed by such a way 
are independent of momentum y (indicated by the horizontal lines which appear in these 
figures) tells the realization of the equilibrium distribution with the temperature T and the 
degeneracy parameters £ v (p). 

In addition to f v and fp are described by thermal equilibrium distribution with the same 
temperature as the photons, we see that the sum of their degeneracy parameters is accurately 
zero: + = (to be more precise, the sum never exceeds 10~ 3 ). This relation is supposed 
to ensure the entropy conservation as discussed above so we should check whether the final 
value of the degeneracy parameter calculated by the numerical simulation is also obtained 
from the conservation laws. The relativistic degree of freedom of fermion with mass m 
(ignoring spin and anti-particle) is 

45 s 45 f, x 2 ( x 2 \ , . 

9s,mass=- 2T - 3 =-^J o ^ exp(e) + 1 (^+3iJ> ( 16 ) 
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FIG. 2: This figure shows the evolution of £p with similar features to Fig^ We see v also keeps 
the equilibrium distribution. Together with Fig. ^ they indicate the relation £„ + £p = holds. 
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FIG. 3: The evolution of the neutrino degeneracy parameter. We plot at y = 7 but it has 
practically no dependence on y as shown in Fig. ^ £p evolves with opposite sign but the same 
absolute value. 
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where e = \J x 2 + (m/T) 2 . Then g s ,muon decreases from 3.11 at T = 100 MeV to 
0.00782 at T = 10 MeV. As for degenerate neutrinos, it is well approximated as g SiV = 
(7/4){(15/7)(£/tt) 2 + 1} and g S)P = (45/2tt 4 )(4 + £)e~ 5 so the latter can be neglected in 
the present case. The other particle species include photon, electron and two types of non- 
degenerate neutrinos (g s ,others = 9). Requiring rji = m/s oc (£ 3 + 7r 2 £) / '{g Symuon + flv(£) + 
g s ,others) to conserve, £ is found to decrease from 10 to 9.526. This reproduces the numerical 
simulation results very well. 

IV. CONCLUSION 

In summary, under the existence of neutrino degeneracy in the early universe, we show the 
total entropy indeed conserves. This is not trivial after the neutrino annihilation process 
which establishes the chemical equilibrium freezes out. We argue that it does as long as 
the neutrinos are kinematically coupled to the rest of the plasma so that the relation + 
£ p = holds. To confirm this argument, we simulate the degenerate neutrino spectrum 
evolution. As the result, we find neutrinos and antineutrinos have the thermal equilibrium 
distributions with the same temperature as the rest of the plasma and with the degeneracy 
parameters satisfying £„ + £ p = 0. Consistently, we see the simulated evolution of the 
degeneracy parameters can be reproduced using the entropy conservation and the lepton 
number conservation. 

Neutrinos with so strong degeneracy that annihilation freezes out before muons annihilate 
have been thought to necessiate some special thermodynamic treatments since Ref. 0] and to 
these days. However, as we discussed in this paper, that is not necessary and the degeneracy 
parameter evolution is calculated using the entropy conservation and the lepton number 
conservation, no matter when their annihilation freezes out. 
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